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Abstract 

We calculate the damping rate for transverse gluons with finite soft mo- 
mentum to leading order in perturbative hot QCD. The internal momenta of the 
one-loop contributing diagrams are soft. This means we have to use effective 
vertices and propagators which incorporate the so-called hard thermal loops. We 
expand the damping rate in powers of the incoming momentum and argue that 
the series ought to converge within a finite radius of convergence. We contrast 
such a behavior with the one obtained from a previous calculation that produced 
a logarithmic behavior, a calculation based on letting the gluon momentum come 
from the hard limit down towards the interior of the soft region. This difference 
in behavior may point to interesting physics around some 'critical' region. 
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1. Introduction 

The quark-gluon plasma (qgp) is a phase of hot hadronic matter that we hope to see in very near- 
future experiments like RHIC and/or LHC. One paramount importance of this plasma is that it allows us 
to see quarks and gluons, if not completely free, at least in a deconfined plasma phase, and get clues on the 
mechanism of confinement. In this respect, many properties of QCD at (high) finite temperature T have 
been investigated [1]. 

One important quantity that is related to the stability of the qgp is the damping rate j(p) of gluons* 
with momentum p in the presumed plasma. This quantity has raised a great deal of controversy in the 
past because, when calculated at p = in a standard loop expansion at finite T, it is plagued with gauge 
dependence both in magnitude and sign [2]. This problem has been solved when recognizing that in such 
conditions, the loop expansion is not necessarily an expansion in powers of the QCD coupling constant g, 
and hence one has to reorganize the perturbative expansion so that one takes into account the resummation 
of the so-called hard thermal loops, which are loop diagrams with hard (i.e., ~ T) internal momenta [3]. 
One then argues that 7(0) is finite, gauge-independent and positive. 

However, a previous calculation [4] suggests that when we let p run from the hard limit down towards 
the interior of the soft region, the boundary of which being physically set by m g , the inverse thermal gluonic 
correlation length or thermal gluonic mass for short (which is of order gT), the damping rate 7t(p) for 
transverse gluons gets a logarithmic behavior ln(l/g). This result is in contrast with the fact that 7 t (p = 0) 
is finite. It is also in contrast with the expectation that the qgp ought to be stable for at least very small 
but nonzero momenta. A similar interest was raised in [5] and a discussion in the context of scalar QED was 
carried. 

In this work, we undertake the expansion in the soft region p < m g of the damping rate "ft(p) in powers 
of p/rrig, i.e., we write: 



in [2] . Our primary aim is to argue that such an expansion is valid within a finite radius of convergence that 
we denote by fi. We do this by explicitly calculating the second coefficient an in the above expansion and 
suggesting that the other coefficients may be calculated in a similar manner. If indeed this expansion is a 
valid one, it would mean that the analytic behavior of the damping rate changes when we cross from the 
region below fj, to the region above. Such a result would suggest that there could be interesting physics to 
investigate in the 'critical' region around fi. 

This paper is organized as follows. After this introduction, we set up the stage in the next section 
for the calculation of the transverse-gluon damping rate. In the third section, we carry the calculation of 

* We discuss electric gluons only. The interesting magnetic sector is more intricate as it is well known. 




(1.1) 



where N c is the number of colors. The quantity 



g 2 N c T 



a t o is just 7t(0) with a t0 = +6.63538 . . ., determined 
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the effective self-enegy to order p 2 , which is the essential new quantity entering the definition of 7t(p), see 
equation (2.9) below. We discuss our results in the last section and finish with few concluding remarks. 



2. Preliminaries 

In this section, we prepare the ground for the calculation of the transverse gluonic damping rate j t (p) . 
We work in the imaginary-time formalism in which the euclidean momentum of the gluon is P p = (po>p) 
such that P 2 = (po) 2 + p 2 where p = p p and po — litnT where n is an integer. After we perform the 
intermediary steps, we obtain the real-time amplitudes via the analytic continuation po = —iu + + where 
u) is the energy of the gluon. The convention we adopt is that a momentum is said to be soft if both u and 
p are of order gT; it is said to be hard if one is or both are of order T [3]. 

We carry our calculation in the Coulomb gauge in which the complete inverse gluon propagator is given 

by: 

D Zu( p ) = p 2 ^ - plxpv - n^(P) + —s^s^pY , (2.1) 

where IT A " y (P) is the gluon self-energy and the last term is due to Coulomb-gauge fixing. In fact, it is most 
suitable to work in the strict Coulomb gauge £c = 0. The gluon self-energy can be decomposed into: 

W V {P) = Sn^(P) + *H^(P) , (2.2) 

where SH is the hard thermal loop and *II is the effective self-energy. P being soft, the hard thermal loop is 
of the same order of magnitude as the inverse free propagator, i.e., <5IT ~ (gT) 2 , while the effective self-energy 
is of an order of magnitude higher, i.e., *n ~ g(gT) 2 . Since the momentum running inside *II is soft, we 
have to use effective vertices and propagators instead of their bare* counterparts when calculating it. This 
ensures the correct expression for the <?(<7T) 2 -correction to the inverse gluon propagator and, in particular, 
that this correction is independent of the gauge. 

The hard thermal loop 5Ti can already be found in the literature, see for example [2,3]. It is real and 
contributes to the determination of the spectrum of the soft gluonic exitations to leading order gT. More 
explicitly, we know that 511 is gauge-invariant and satisfies the identity P^5Yl^ u (P) = 0. This means its 
components can be expressed in terms of only two independent scalar functions denoted by 511; (P) and 
SU t (P) such that: 



SU m (P) = SU t (P) ; SU 0t (P) = -?°P 6^ (P) ; 

P 

sn^(p) = (s ij -ff) sn t (P) +fp> m.i{P) 



(2.3) 



The expressions of 5Hi (P) and <5IT t (P) read [3] : 



* 'Bare' refers here to the usual quantities one considers as dictated by the Feynman rules. 
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where the Q n are Legendre functions of the second kind. As already mentioned, m g is the gluon thermal 



mass and, to lowest order, is equal to ^JN C + (1/2) Nf gT/3, where Nf is the number of flavors. 

The effective propagator for soft gluons that intervenes in the calculation of the effective self-energy is 
obtained by inverting (2.1) while disregarding *II. In the strict Coulomb gauge, its nonzero components are 
*A£?(P) = *Ai(P) and *Ag(P) = (<J« - p l p>) *A t (P), where *A t and *A ; are given by: 

* A ^= p 2 -W * Ai(p) =^k(p)- (2 - 5) 

After analytic continuation to real energies, the pole in uo of *A t ^ yields the dispersion relation uj t (i)(p) for 
the transverse (longitudinal) gluons to order gT. One finds for the soft transverse ones: 

[,3/P\ 2 9 / p \* 704 / p n6 91617 / p \« i 

As mentioned earlier, the hard thermal loop 511 is real above the light cone, and so the poles of (2.5) 
are real, which means the gluons are not damped to this order gT; this is clear from (2.6). In order to 
get the leading order of the damping rates, we have to include in the dispersion equations the contribution 
from the effective self-energy. *II has of course a more complicated structure than 5U. It satisfies the less 
restrictive identity P^ "IP" \P) P v = [3]. This means that in general, ' (P) can be written in terms 
of three independent scalar functions, but in the strict Coulomb gauge, only two of these are relevant. The 
transverse dispersion relation including the self energy reads: 

-ft? + P 2 - m t (-in t , P ) - *u t (-m t ,p) = o , (2.7) 

where *II t is given by: 

*n t (p) = - fp>) *u l \p) . (2.8) 

The transverse gluon damping rate is defined by 74 (p) = — Im ftt(p). Since it is 5-times smaller than the 
energy totip), we can write from (2.7): 

f ) = Im * n t(~^ ; -P) / 2 

lt{P ' 2cv + £SIl t (-iuj,p) u=u t ( P )+M+' 
The denominator in (2.9) is easy to get since we already have an expression for the hard thermal loop <5II t . 
Indeed, we have: 

2w t ( P ) + d u Sn t (-iu J ,p)\ u=Ut(p)+i0+ =2[l + ^(-^) 2 + ...] . (2.10) 
This means our main task is to calculate the imaginary part of *Tlt- This we do in the sequel. 

3. The imaginary part of the transverse effective self-energy 

In the Coulomb gauge, the only diagrams that contribute to the imaginary part of the effective self-energy 
above the light cone are the three-gluon and four-gluon one loop-diagrams with soft internal momentum [3] . 
Hence we write: 

2 AT" 

Im *lT" y (P) = - L^. Im Tr soft [*r^ A<T (P, -P,K, -K) *A$?(K) 

+ *T a ^{-Q,P,-K) *A^'(K) *r x '^'(-K,P,-Q) *A£>(Q)], (3.1) 
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where K is the internal loop-momentum, Q = P — K and Tr = TJ2 ka J jk^s • The subscript 'soft' means 
that only soft values of k are allowed in the integral. Eq (3.1) is what one would normally write for the three- 
gluon and four-gluon contributions to the imaginary part of the gluon self-energy, except that everywhere, 
bare quantities are replaced by the corresponding effective ones. 

We have already given in the last section the expressions of the nonvanishing components of the effective 
propagators (see eq (2.5) and the text before it). The gluon effective vertices can be written as: 

* r («) = r (n) + §T (n) . n = 3j4j (3.2) 

where the first term is the QCD gluon tree vertex and the second one sums up the contributions from hard 
thermal loops with n external legs. In the case n = 3 it can be written as: 

^ ( _ g)Pi _^ ) = _^A ( _^ p> _ 0)=3m ,^^(^_^) t (3 . 3) 

where S = (i, s) and fls is the solid angle of the unit vector s. Also, PS = ipo + ps, etc. In the case n = 4 
we have: 

xr iw\*,p p v v\ o 2 f d ^s S fl S"S x S a ( ipp - ik ipo + ik \ 

ST (f , -P, K, -K) = 3m g J — pSRS (^^-^--^-^j. (3.4) 



(3.5) 



To be complete, we give the expression of the three-gluon tree vertex: 
T^ X {-Q,P,-K) = -T X ^(-K,P,-Q) 

= (P + KY S vX + {Q- K) v 5 X » -(P + Q) x 5^ , 
and that of the four-gluon tree vertex: 

V nv\o -P^k,-K) = 2 8^5 Xa - 6^ X 6 V " - 5^5 uX . (3.6) 

From eqs (2.8) and (3.1) above, we can write more explicitly the expression of the transverse effective 
self-energy: 

Im *II t (P) = (5^-Pf)lmTj2 I 7^ [*r^ 00 (P, -P, K, -K) *A t (K ) 

+ *r ijm "(P, -P, K, -K) (S mn - k m k n ) *A t (K) 

+ *T 0M (-Q, P, -K) *A l {K) *T 0j0 (-K, P, -Q) * A Z (Q) (3 7) 

+ *T Mm (-Q, P, -K) (S mn - k m k n ) *A t (K) *T nj0 (-K, P, -Q) *A l {Q) 

+ *T mM (-Q, P, -K) *A l {K) *T 0jn (-K, P, -Q) (S nm - q n q m ) *A t (Q) 

+ *T mir (-Q, P, -K) (S rs - k r k s ) *A t {K) *T sjn (-K, P, -Q) (6 nrn - q n q m ) *A t (Q) ] . 

There are six contributions: two from the four-gluon vertices (subscript Ag in the sequel) and four from the 
three-gluon vertices (subscript 3g). Each one has to be calculated separately. As an illustration, we show 
how we carry the calculation corresponding to the contribution from the three-gluon vertices where the two 



effective propagators involved are both longitudinal. We denote this contribution by lm*U t 3 g ii(P) and the 
others correspondingly. The other contributions are manipulated in a similar way, with varying difficulties 
that we comment on later in this section. Also, we will take henceforth m g = 1 and all momenta and energies 
are in units of it. This simplifies considerably the expressions we write down and, if and when needed, the 
m g — dependence can be recovered in the final results. 

Using the expressions of the gluon vertices we gave in eqs (3.2)-(3.6), we write: 



Im*n t3s H(P) 



9 2 N C 
' 8tt 2 

+ 6 

+ 9 



fc 



d 3 k 



ImT E/^ (p- 2k ) -(p- 2k )p : 



dtls (p — 2k)s (p 2k)p sp / ik iq \ 
PS \~KS ~ QS) 

dQ S2 s i s 2 -sips 2 p / iko 



4tt PS! PS 2 



( ikg _ m_\ ( 
\KS X QSJ\ 



ik 



iqo 



KS 2 QS 2 , 



t A l (K)*A l (Q) . 

(3.8) 

We first work on the term that does not involve a solid-angle integral. We use the relation (p — 2k) 2 — (p — 
2k)p 2 = 4/c 2 sin 2 ip, where tp = (p, k), and integrate over the solid angle of Kk after we expand the effective 
propagator at the momentum Q = P — K in the following manner: 



*A,(Q)= [l- P cos^ fe + ^(^^d fc +cos 2 ^ 2 ) +...] *A l (q ,k) 



(3.9) 



where dk = d/dk. We find: 

<f'k 

feo 



ImT E / ^[(p-2k) 2 -( P -2k) P 2 ]*A i (X)*A i (g) 



(3.10) 



*A l (q ,k) 



The next step in to perform the sum over fc , but we take care of that a little later. 

We leave the expression in (3.10) as it is for the moment and turn to the term in (3.8) that involves one 
solid-angle integral. It is sufficient to concentrate only on the piece that contains the ratio iko/KS because 
the other one that contains the ratio iqo/QS is in fact equal to the first one. We need an expression for the 
solid-angle integral. For this purpose, it is best to mesure the solid angle Q$ = (0, 4>) with respect to k such 
that 9 = (k, s). Also, we expand l/PS in the following manner: 



PS inn V Wn Bn / 



(3.11) 



Wo Po 

an expansion valid as long as \p/ipo\ < 1, which is satisfied for soft gluons before and after analytic contin- 
uation. With this, we can write: 



/ 



dQs (p — 2k)s — (p — 2k)psp ikn 

~KS 



4tt 



PS 

-2ksmip — - 
Wo 



d£ls sin ip cos 9 + cos tp sin 9 sin 



47T 



(iko + kcos9) 
6 



V ipo pi ' 



(3.12) 



With the relation ps = cos ip cos 9 — sin ip sin 9 sin </>, the angular integrals are performed straightforwardly 
and we end up with: 



/ 



47T 



PS 



dfl s (p - 2k)s - (p - 2k)p sp ik Q ik 2 ^ fji ik o ^ \ , P ( Jk Q 
tt^ = - - — (1 - x ) 2 1 - -r- Qok + — x 3— 



^l-2, 2 ) + (l-5* 2 )| 



(3.13) 



+ ... 



where a; = cos?/> and Qofc = Qo{-jf)- We plug this expression back into the integral over d 3 k and use (3.9) 
to perform the integral over the solid angle of k. That is quite straightforward and we get: 



6 ImT 



^ f d 3 k f dn s (p - 2k)s - (p - 2k)psp / ik iq Q \ ,^* A {n , 

l^J^J^- ps [ks-qs) mk) AliQ) 



k 

= -16 ImT 



r c 



k 2 dk*A l (K)^\(l 



ik 



Q()k 



Spl 



k 
+ 



ink 



(3.14) 



+ ? ( 3 x + i 1 + 3 5l) Qok ) dk +p »{ 1 - x Qofe ) (l dk + \ dl )\ + ■ ■ ■ } * Al{qo > k) ■ 

Here also we leave this expression as it is for the moment and turn our attention to the term in (3.8) 
involving two solid-angle integrals, this term is actually equal to 



feo 



dilsi f d£ls2 sis 2 - sips 2 p 



4tt PS1PS2 \KS 1 KS 2 KSxQS 



0Q0 \ 
1QS2) 



and we first work out the piece that contains k 2 /KS\ KS 2 - In order to be able to carry forward, it is best 
to write the double solid-angle integral in the following manner: 



dfl Sl 
An 



d£ls 2 Si§ 2 - §ipS 2 p fcg 2 

4tt PSi PS 2 KSi KS 2 ~ 



sin 2 



dfls cos 9 
~IV PSKS 



+ cos ip 



dfla sin 9 sin < 



4tt PSKS 



+ 2 cos V> sin 



dQ.s cos 9 f dCts sin# sm< 



, f dCls cos# f 
^ J ^ PSKS J 



4tt PSKS J 4tt PSKS 



(3.15) 

Each single solid-angle integral that is involved in the above expression can be worked out straightforwardly 
as before, using the expansion (3.11). Putting things together, we find: 



d£l S2 §iS 2 - S!pS 2 p 



kl{\ - X 2 ) 



4tt PSi PS 2 KSi KS 2 

px ( ^ ik „ \ / _iki 
Wo 

ik „ \ f 2 4 



Plk 2 



^+('+i)«»y 



+ .. 



(3.16) 
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We put this expression back under J d 3 k and perform the integral over the solid angle of k. We get: 



k 2 



3 tr * * 



dtts 2 §is 2 - §ips 2 p 
4tt PSi PS 2 KSi KS 2 

\ 2 P 2 fl (iko 



x ( s *. + ( l+s |) <3h )*.ffi( I . W (^ + « 



2 /4 



*Ai(q ,k) . 

(3.17) 

Next we turn to the piece that contains fcoQo / QS 2 ■ To get a manageable expression for the double 
solid-angle integral, it is most suitable to measure the solid angles with respect to p. Then we have: 



dQ 



Si 



A.-K 



dn S2 §is 2 - §ips 2 p k q 



4tt PSi PS 2 KS 1 QS 2 



koqo 



dfls sin 9q sin c/)q f dfls sin 9q sin </>o 



4tt PSKS 



4tt PSQS 



(3.18) 



where (6>o,0o) is the solid angle of s with respect to p. Each solid-angle integral can be carried through 
separately. We get the following result: 



/ 



cifis sin O sin </> sin ip 
PSKS 



>->^( 3 X + ( 1 + ^ 



ip k 
P 

IPI 



1 r l 
3 



(1 - 2x 2 ) + (1 - hx 2 )^ (l - 3x 2 + (1 - 5x 2 )^) ^Q ok 



(3.19) 



and a similar one for the other integral in (3.18) where one replaces K by Q and the angle (p, k) by (p, q). 
We then multiply the two obtained expressions and put back the result under T^2 kg J d z k/A-K. But here 
the integral over the solid angle of k is not straightforward yet because the integrand still depends on q and 
the angle (p,q). Hence a further expansion is necessary, but instead of expanding Qo q directly it is most 
suitable to write it in terms of *A ; _1 (Q) using (2.4) and (2.5) and expand the resulting expression. The 
calculation carries thereon straightforwardly and we get: 



18ImT^ 



d 3 k, 
47T 



A ; (tf)*A ; (Q) 



*A,(K) 



-Po 



ikoiqok 2 



dnsj_ f dn S2 Sis 2 -Sips 2 p fc <7o 
t J 

rH r 1 



4tt P5iPS* 2 KS 1 QS 2 3 



= - ImT E 



k z dk 



5P 2 



2ik iq k 2 + ^(3 + k 2 + 3k 2 ) (3 + k 2 + 3qg) 



+ ... 



- iq ip k(3 + k 2 + 3kl)d k + ik a iq p 2 Q k 2 (^d k + ^k) 

(3.20) 

It is time now we get into performing the sum over kg. Because of the complicated fco— dependence of 
the expressions we work with, it is best to perform the sum using the spectral representation of the different 
quantities involved [6] . For example, we have for the effective propagators: 



dre lkoT I dup t ^(u,k)[l + n(u)]e-^ T , 

J —CO 



(3.21) 



where n(ui) = l/(exp(u;/T) — 1) is the Bose-Einstein distribution and the spectral function pi it is given by: 
Pt,i(w, k) = [S(oj - Lu t ,(k)) - S(uj + wt,i(fc))] + fc) 6(fc 2 - w 2 ) , (3.22) 

expression in which the residue Z t j(k) is given by: 

zm - ^ 2 -fc 2 ) , i 4^) f3 23) 



and the cut f3 t .i reads: 

?M{k 2 - u 2 ) 



4fc 3[ (fc 2 _ QJ 2 + + fei_^ ln |±H ))2 + ( 3™ (fc2 _ w2)) 2] > 

3lo 



(3.24) 



2fc[(3 + fc 2 "tln^) 2 + (^) 2 ] ' 



After we replace the effective propagators and other similar quantities by their spectral representations, we 
perform the integrals over the imaginary times. One such integration yields a delta-function and the others 
energy denominators. Actually we arrange our expressions in such a way that we have always only two 
imaginary-time integrals, which ensures that we get only one energy denominator. When this is performed, 
we analytically continue ipo to the on-shell real energy u>t(p) + *0 + . The imaginary part of the different 
contributions is obtained using the well known relation l/(x + i0 + ) = Prl/x — inS(x). We apply this 
technique to the expressions in (3.10), (3.14), (3.17) and (3.20) and sum up the results to obtain: 



T *rr ,™ 9 N cT f°° „ f j f j 6w 2 p 2 

Im U t3 gii(P)=— / dk / diox / dw 2 P11P12 + -, — ©l P12 + -=■ 

24tt J q J_ 00 J_ 00 Lluilu 2 kuj 2 5 



3fc 



2 



x (3 + U* - l^ )pa Pa + £ (l - ||) (l - 9 |i) 9, „ + (13 + 3* 



+ lOwi - 9cjI) pa d kP i 2 + 6— (l + 2^ - 3^-) 6i d kPt2 + (2 + 5wi) pn 0£p I2 

3w 2 2 54fc 4 18^ 2 i i , , 

H rBi4/? i2 Piipi20 ui J-Oipno^ +... 5(1 -wi -w 2 J , 

W 2 K WiW 2 W 2 K J J 

(3.25) 

where p/j denotes pi(u>i, k) with i = 1,2, and 9i = 6(fc 2 — w 2 ). In (3.25), we have used n(wj) ~ T/wj since 
only soft values of u>i are to contribute to the integrals. Also, we note that many terms drop out in the end 
because they do not hold an imaginary part. The remaining integrals in (3.25) are to be carried through 
numerically, more on this in the last section. 

As we said earlier, the other contributions to Im*IIt(P) in (3.7) are manipulated along the same lines. 
Concerning the 3g-contributions, the intermediary steps are much longer, and at times quite intricate. For 
example, one thing that doesn't look obvious from the outset is that, when performing the sum over ko, it 
is necessary to group the different terms in a certain manner such that expressions containing fco in a way 
that prevents us from going to the spectral representation do cancel out and the technique can then be used. 
Another grouping of terms is also necessary in order to ensure that we end up at each time with only one 



energy denominator. Such subtleties do not occur in the 11 contribution we worked out above. It would be 

too long and quite cumbersome to report on the intermediary steps of these contributions, and so we content 

ourselves with giving the final resutls, similar to (3.25). For the It and tl contributions, we find: 

t *tt ( r>\ t *tt rm 9 2 N C T1 f° „ r +co cLj 1 f +QC duj 2 
Im U t3gU (P) =Im H t3gH (P) = 



pOO /"+oo 

dk 

JO J -co 



duii 



3wi 



24tt 2 

(k 2 -u> 2 fQ 1 p t2 + ^(k' • 9 i /*/.. + 



oo w 2 
^2 



-i&k 2 (k 2 -wiypa pn 



36 



20 - 12A- 2 + 93A- 4 + 50fc 6 
/72 



+ ( - ^ + g + 12fc 2 + 78fc 4 ) Wl + (1= - 80 - 2fc 2 - 186fc 4 ) w 2 + g + 36 - 176fc 2 ) u; 3 

+ ( 



84 
fc 2 



3 2 



- 179 + 222fc 2 
15\ 



Ptl PZ2 + 



12 



©lPi2 



- 30/ccJi 



+ — w? — — w? 

fc A: 3 



©i p;2 + [-12k - 72k 5 + (-24k - 12k 3 + 48fc 5 )^i + (-12k + 144fc 3 )w 2 

+ (12k - 96k 3 )ujl - 12kuj\ + 48fcwf ] p tl d k p l2 + [69k + 14k 3 - 3k 5 + (-54k + 8k 3 + 6k 5 )ui 

+ (2k + 2k 3 )ul + (8fc - l2k 3 )uo 3 - \\kuj\ + Qkuo\ ] p n d k p t2 
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9 „ 9 3 , 
-wi + 9^ 2 + ^w? 



18 4 3 5 9 6 
- ¥^ + + k^ 



3 18 4 



©1 <9 fc p; 2 



©i <9 fc p t2 + [9wi - 18w 2 
+ [-8fc 2 - 16fc 6 - 16fcV + (-8fc 2 + 32A- 4 )u 2 - 16A- 2 w 4 ] p tl d\p l2 + [30fc 2 + 4fc 4 - 2fc 6 
+ (-24fc 2 - 8/c 4 Vx + {-4k 2 + 4k A )ui\ + 8k 2 us 3 - 2k 2 ui\ ] p n d 2 k p t2 



3^i(fc 2 - c 2 ) 2 S 1 d 2 kPt2 + 6^(fc 2 - u() 9! d 2 p ;2 + 54fc 2 (fc 2 - ul) 2 p n Pl2 a, 



A: 3 



+ 9^| (fc 2 - c 2 ) 2 9! p t2 d Ul - 18^ (fc 2 - c 2 ) 0! Pl2 d ui 



(5(1 - ui — uj 2 ) 



(3.26) 

The above expression is rather long, partly because these two contributions do not benefit from a possible 

symmetry between u>i and uj 2 , which can lead to some useful simplifications. Indeed, though the algebra is 

by far the most tedious for the tt contribution, the use of this symmetry renders its final result relatively 

simpler. It reads: 

g 2 N c T 



Im*U t3gtt (P) =- 



24tt 

4 

+ ~4k 



dk 



+oc dLJ 1 

-oo 



+O °dco 2 



LU 2 



36(fc 2 + ujiuj 2 ) 2 pn p t2 - 6 "p-( fc2 - w i) ©i Pti 



1 



ISA; 2 + 122A 4 + ±(-728 + 854fc 2 ) unu 2 + (J£L - ™)lo 2 lo 2 - ^ 3 co 3 



Ptl Pt2 



- 5fc 2 - 16wi + 2w 2 
'48 



21 4 



©lPt2 

108 



32k 3 + 20k b 



(112k + 124/c 3 ) W! + - 36/c - 116/c 3 )w 2 - + 232/c)w 3 + + 156fc)w 4 



108 5 

H w? 

k 1 



60 



Pti a fc p t2 - (l - (k 2 + 3CO! - 3co 2 ) 0i d kPt2 



+ [8/c 4 + 20fc 2 (l + A: 2 ) ui + 12(1 + A- 2 ) - (12 + 32k 2 ) uj\ - \2uj\ + 12ojf } p t i d 2 k p t2 

UJ 3 

- 3-rj(fc 2 - uj 2 ) 9i d 2 kPt2 - 108(A- 2 + uilo 2 ) 2 p tl p t2 d Ul 



+ 18^i(A; 2 -^ 2 )9ip t2 ^ 1 



+ ... 



(5(1 - wi — w 2 ) . 



(3.27) 
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The two 4g-contributions are easier to work out and can be cast in different ways depending on how we 
perform the expansion. One compact way is to write: 



■2]\tt< r°° r+°° r+°° u afu^ _ ,.,2 

2 



a 2 N t f°° r +aa r +oc 

Im*n t4gl (P)+lm*n t4gt {P) =^4^ / dk J dwi / ^ 

24fc 3 . ilnWl . , t , 6(4^-1) 



„ k n 6(k 2 -uf) „ 
6— 01 P 12 + —. - 01 pt2 

W2 KW 2 



5i p/2 + 12 — e(wi) (5i p t2 H : 01 p t 2 

WIUJ2 Ul2 KUJ2 



(3.28) 



6/cwi r . 18A; _ 
H e(wi) d Wl [<)(« - wi) - + wijj H 01 P12 o u 

U>2 1^2 

_18(^ -| + _ ] ^_ 

KW2 - 1 - 1 

where 5i denotes <5(fc 2 — w 2 ) and e(u!\) is the sign function. 
4. Discussion 

Recall that our aim in this work is to calculate to order p 2 the damping rate for soft transverse gluons in 
QCD at high temparature T. We have mentioned in the introduction that our main motivation is to suggest 
that the analytic behavior of this damping rate, and perhaps similar quantities, may change around some 
scale p, see the first section. 

The damping rate 7t(p) is given by equation (2.9). The denominator is given to order p 2 in (2.10). The 
numerator Im*n t (P) is the sum of the expressions given in (3.25)-(3.28). As we said in the last section, its 
actual value is to be determined numerically, an issue we comment on a little later. But at least from an 
analytic point of view, our calculation shows that the expansion in powers of p 2 of the damping rate is quite 
feasible: we recover the order zero in p 2 given in [2] and find an expression for the next one. Furthermore, 
our steps show that the calculation of the higher orders proceeds in a straightforward, though practically a 
lot more tedious, manner. 

Supposing for a moment that the numerics go smoothly, it is interesting to ask about the nature of this 
scale p. The first suggestion that comes to mind is that p is related to magnetic effects which are believed 
to manifest themselves at the next natural order beyond m g , i.e., g 2 T. If this is so, there are two interesting 
points to mention in this regard. First of all, it would be interesting to understand how the magnetic effects 
get into play with respect to other effects and how and why they make this damping rate and possibly other 
similar quantities change their analytic behavior. Also, in the case of this change being indicative of some 
sort of phase transition, it would be interesting to understand the nature of this transition and the 'critical' 
behavior of QCD around it. 

However, it may be that the scale p is due to some other effects. Though this maybe remote a possibility, 
it would certainly be very interesting to contemplate into the matter. One way to start looking into this is 
to notice that the coefficients au in the expansion (1.1) are all pure numbers. Thus, with more coefficients 
calculated and assuming that they are all finite, it is possible to determine approximately the radius of 
convergence of the series, using for exemple a Pade approximant technique, and hence determining p with 
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respect to m g . If by other independent means we have an expression for the running coupling constant g, 
we can compare the approximate value of n/m g to this latter and carry a discussion thereon. 

But at this stage, all this is still speculative. Indeed, though, as we said, the analytic expression of the 
effective self-energy to order p 2 is quite clean, it is not apriori obvious it is so as far as its numerical value 
is concerned, it is true that the coefficient a t o we get from our calculation has a finite value, but this is not 
certain for a t \- First of all, we have expressions involving products of derivatives of delta- functions we have 
to clarify precisely how to carry through with. But more importantly we think is that we are not assured 
of an infrared safeness with regard to the integration over the internal soft momentum k. This issue needs 
careful and thorough investigation, which is beyond the scope of the present work. 

From a physical standpoint, the fact that ato is finite means that the quark- gluon plasma is stable 
against zero-momentum gluonic excitations. For practical purposes, it is sufficient to take 7t(p) ~ 9 2 ^ T ato 
for small enough gluonic momenta. This implicitly presupposes that for such small non-zero momenta, the 
plasma remains stable. But this means that, QCD being the adequate theoretical framework to describe the 
qgp, 7t(p) ought to have a smooth analytic behavior for such small momenta, and hence the highcr-than-zero 
coefficients in our expansion (1.1) should be expected to be finite. This is the basic motivation of our work. 
If on the contrary they happen to be not so, because of infrared divergences for instance, then one must look 
into the matter more closely and try to reconcile between what we expect and what we get. 

After these issues are settled should come the discussion of the physical implications of these results. 
For example, the sign of ato being positive, is the sign of a t \ also positive? If not, then at what momentum 
instabilities start to apprear? We think all these issues are worth pursuing. 
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